AFRL-VA-WP-TP-2002-304 

CONCEPTS  FOR  CONSTRAINED 
CONTROL  ALLOCATION  OF 
MIXED  QUADRATIC  AND  LINEAR 
EFFECTORS 

David  B.  Doman 
Andrew  G.  Sparks 


APRIL  2002 


Approved  for  public  release;  distribution  is  unlimited. 


This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright 
protection  in  the  United  States. 


20020912  082 

AIR  VEHICLES  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7542 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
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Abstract 

Concepts  for  new  constrained  control  allocation  strate¬ 
gies  are  developed  that  deal  with  systems  where  mo¬ 
ments  are  nonlinearly  related  to  effector  deflections 
such  as  those  encountered  in  the  case  of  yawing  mo¬ 
ment  contributions  from  left-right  effector  pairs  on 
aircraft.  These  concepts  are  illustrated  by  consider¬ 
ing  single  and  multiple  left-right  pair  effector  mixing 
problems  for  moments  that  lie  in  the  roll-yaw  moment 
plane.  Methods  for  generating  the  boundary  of  an  at¬ 
tainable  moment  set  for  a  class  of  multiple  non-linear 
effectors  and  for  clipping  unattainable  moment  com¬ 
mands  with  axis  prioritization  are  presented. 


1  Introduction 

Left-right  aerodynamic  surfaces  such  as  ailerons,  flaps 
and  elevons  on  aerospace  vehicles  often  have  non¬ 
linear  contributions  to  the  total  vehicle  yawing  mo- 
nfient.  Since  these  surfaces  are  not  normally  used  as 
primary  yaw  axis  effectors,  these  non-linear  contribu¬ 
tions  are  often  ignored  and  their  effects  are  treated  as 
disturbances  that  are  rejected  by  a  feedback  control 
system.  In  cases  where  the  effectiveness  of  the  pri¬ 
mary  yaw  axis  effector  is  degraded  due  to  failures  or 
as  a  result  of  being  in  the  wake  of  the  body  at  high 
angles  of  attack,  these  secondary  nonlinear  effects  be¬ 
come  important. 

Numerous  control  allocation  and  control  effector  mix¬ 
ing  algorithms  have  been  developed  over  the  past 
decade  and  excellent  survey  papers  have  been  writ¬ 
ten  that  point  out  the  strengths  and  weaknesses  of 
the  existing  approaches[l,  2].  These  control  mixing 
and  control  allocation  algorithms  are  capable  of  deal¬ 
ing  with  systems  where  the  moments  were  linearly  re¬ 
lated  to  control  effector  positions  and  had  the  ability 
to  account  for  constraints  on  those  positions.  Some 
of  these  algorithms  generate  constrained  control  effec¬ 


tor  commands  that  ensure  that  the  effectors  are  never 
driven  beyond  their  physical  limits.  Most  of  the  algo¬ 
rithms,  however,  assume  that  linear  relationships  exist 
between  the  pseudo-commands  (i.e.  controlled  vari¬ 
able  commands)  and  the  effector  positions.  While  this 
assumption  is  at  least  locally  valid  for  many  of  the  con¬ 
trol  surfaces  on  found  on  aircraft,  there  are  exceptions. 

One  particular  case  where  this  assumption  can  result 
in  incorrect  control  surface  deflections  and  return  un¬ 
necessarily  conservative  results  involves  the  use  of  left- 
right  aerodynamic  control  surfaces  on  aircraft.  Ex¬ 
amples  of  left-right  aerodynamic  surfaces  include  left- 
right  elevators  and  left-right  ailerons.  This  type  of  sur¬ 
face  can  generally  produce  pitching,  rolling  and  yaw¬ 
ing  moments.  While  these  surfaces  normally  produce 
pitching  and  rolling  moments  that  are  locally  linear  in 
control  surface  deflection,  they  can  have  a  highly  non¬ 
linear  contribution  to  the  yawing  moment  especially 
when  parasitic  drag  dominates  induced  drag  effects. 

In  particular,  these  surfaces  can  generate  yawing  mo¬ 
ments  that  are  of  the  same  sign  whether  they  are  de¬ 
flected  up  or  down.  This  is  because  they  generate  a 
drag  force  on  the  side  of  the  vehicle  on  which  they 
are  located  regardless  of  whether  they  deflected  in  the 
positive  or  negative  direction.  The  yawing  moments 
generated  by  these  effectors  are  usually  small  when 
compared  to  a  primary  yaw-axis  effector  like  a  rudder; 
however,  their  effects  can  become  significant  when  a 
rudder  fails  or  when  the  aircraft  is  operating  at  high 
angles  of  attack  where  flow  over  the  rudder  is  inter¬ 
rupted  by  the  body. 

This  particular  nonlinearity  is  used  as  motivation  and 
serves  as  an  example  for  the  concepts  explored  in  this 
paper.  Methods  for  generating  the  boundary  of  an  at¬ 
tainable  moment  set  for  a  class  of  multiple  non-linear 
effectors  and  for  clipping  unattainable  moment  com¬ 
mands  with  axis  prioritization  are  presented.  The 
techniques  presented  are  applicable  to  more  general 
nonlinear  control  allocation  problems. 
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2  Attainable  Moment  Set  for  a  Single 
Left-Right  Effector  Pair 

Durham[3]  developed  and  subsequently  refined  meth¬ 
ods  for  determining  attainable  moment  sets  (AMS) 
for  effectors  that  generate  moments  M  that  are  linear 
functions  M  =  of  the  effector  positions  6  subject 
to  constraints  on  those  positions  S  =  <  ^  <  ^}‘ 

The  general  solution  to  finding  the  attainable  moment 
set  for  an  over-actuated  linear  system  with  position 
constraints  involves  constructing  a  polyhedron  in  mo¬ 
ment  space.  Potential  vertices  are  are  constructed  by 
locking  all  control  effectors  at  their  extreme  positions 
in  all  possible  combinations  while  allowing  two  effec¬ 
tors  to  traverse  the  range  of  their  possible  positions. 
Durham’s  algorithm  connected  the  vertices  to  form  po¬ 
tential  boundary  facets  and  determined  which  facets 
were  on  the  boundary  of  the  AMS. 

Left-right  (LR)  pairs  of  effectors  such  as  ailerons,  el¬ 
evators  and  flaps  as  mentioned  previously,  generate 
non-linear  contributions  to  the  vehicle  yawing  mo¬ 
ment.  Figure  1  shows  the  yawing  moment  that  is 
generated  by  deflecting  the  right  eleven  of  a  partic¬ 
ular  lifting  body  vehicle.  One  can  see  that  a  parabolic 
fit  provides  an  adequate  approximation  of  the  original 
data.  Generation  of  linear  fits  to  this  data  for  use  in 
a  conventional  control  allocator  that  assume  a  linear 
relationship  between  the  moments  and  control  effector 
deflections  is  problematic.  This  is  because  lines  fitted 
using  either  negative  deflection  or  positive  deflection 
data  results  in  lines  with  slopes  of  opposite  signs.  Fur¬ 
thermore,  since  no  single  line  can  accurately  model  the 
data  and  no  matter  which  line  is  selected,  the  sign  of 
the  yawing  moment  estimate  will  be  incorrect  half  of 
the  time. 


xlO^ 


Figure  1 :  Yawing  Moment  due  to  Right  Elevon  Deflection 


Here  we  examine  a  case  where  a  LR  pair  of  elevens 
exerts  yawing  moment  N  that  is  proportional  to  the 
square  of  deflection  and  rolling  moment  L  that  is  a 
linear  function  of  deflection 

N  =  Nie^Sf^  + 

£j  z=:  Li^Slc  “h  Lre^rC' 

From  the  curve  fits  to  aerodynamic  data  for  a  lifting 
body  model  at  a  subsonic  flight  condition,  we  select 
=  -Nie2  =  34/^  -  Ib/deg^  and  Lre  =  -Ue  = 
-4610/t  “  lb /deg.  The  AMS  is  generated  by  holding 
each  surface  at  a  fixed  deflection  while  allowing  the 
other  surface  to  vary  over  its  range  of  possible  values 
and  is  shown  in  Figure  2. 


Figure  2:  Nonlinear  AMS  and  Lines  of  Constant  Deflec¬ 
tion  in  the  Roll- Yaw  Moment  Plane 

The  AMS  is  also  shown  for  the  yawing  moment  as  a 
linear  function  of  the  deflections.  For  the  linear  ap¬ 
proximation  we  fit  data  over  half  of  the  range  of  deflec¬ 
tions  to  obtain  =  ±676/t  —  Ib/deg.  The 

sign  ambiguity  results  from  the  yawing  moment  deriva¬ 
tives’  dependency  on  the  sign  of  the  deflections;  either 
choice  of  sign  will  only  yield  a  yawing  moment  with 
the  correct  sign  over  half  of  the  deflection  range.  The 
AMS  for  the  linear  case  shows  that  a  control  allocator 
that  incorrectly  assumes  a  linear  relationship  between 
the  yawing  moment  and  the  effector  deflections  will  be 
constrained  to  rolling  and  yawing  moments  that  lie  on 
one  of  two  lines  in  the  roll-yaw  plane.  The  AMS  for 
the  improved  effector  model  shows  that  larger  regions 
in  the  roll-yaw  plane  are  reachable  and  the  intersec¬ 
tion  of  lines  of  constant  effector  deflection  define  the 
proper  blend  of  effectors  required  to  achieve  a  partic¬ 
ular  rolling  and  yawing  moment  that  lies  within  the 
AMS. 
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3  Clipping  Infeasible  Moment  Commands 

It  is  interesting  to  note  that  it  is  not  possible  to  gener¬ 
ate  yawing  moments  without  generating  a  correspond¬ 
ing  rolling  moment  with  one  set  of  left-right  effectors. 
This  point  begs  the  question  of  how  to  clip  rolling 
and  yawing  moment  commands  that  lie  outside  of  the 
AMS.  When  a  moment  command  requires  violation  of 
an  effector  position  constraint,  Durham[4]  suggested 
that  the  direction  of  a  moment  command  should  be 
preserved  and  that  the  command  should  be  clipped 
such  that  the  commanded  moment  vector  touched  the 
edge  of  the  AMS. 

For  the  case  under  consideration,  one  can  see  that  pre¬ 
serving  the  direction  of  a  pure  yawing  moment  com¬ 
mand  will  result  in  the  command  being  clipped  to  zero. 
A  proposed  method  for  clipping  commands  that  lie 
outside  of  the  AMS  that  provides  the  ability  to  pri¬ 
oritize  moments  is  given  below.  When  a  command  is 
infeasible,  we  find  a  way  of  choosing  a  point  on  the 
boundary  of  the  AMS  that  is  closest  in  some  sense  to 
the  infeasible  command. 

We  propose  to  minimize  the  sum  of  weighted  square 
distances  between  the  AMS  boundary  and  the  desired 
point  in  moment  space 

min  [{La  -  Lf  +  w{Na  -  iV*(L))2] ,  (2) 

where  is  a  weighting  factor  that  can  be  used  to 
weight  the  relative  importance  of  achieving  the  rolling 
or  yawing  moment  and  Nb{L)  is  the  equation  for  the 
AMS  boundary  which  is  closest  to  the  desired  yawing 
moment  The  first  order  necessary  condition  for  an 
extremum  of  the  cost  function  is 

L  -  Lrf  +  wNl,{L){Nb{L)  -  Nd)  =  0,  (3) 

which  can  be  solved  for  L  either  numerically  or  an¬ 
alytically  in  special  cases.  Figure  3  shows  the  effect 
of  varying  weights  in  the  cost  function  for  a  hypo¬ 
thetical  infeasible  moment  command.  Note  that  when 
w  ^  0  the  desired  rolling  moment  is  achieved;  as  is 
increased,  the  error  in  yawing  moment  is  reduced  to  as 
small  of  a  value  as  possible  while  the  rolling  moment 
error  increases.  This  type  of  strategy  may  be  useful 
if  other  effectors  are  available  to  reduce  the  undesired 
rolling  moment  if  a  particular  left-right  pair  must  be 
used  as  a  primary  yaw-axis  effector. 

To  account  for  nonlinear  effectors  one  must  consider 
the  most  general  form  of  the  control  allocation  prob¬ 
lem.  This  problem  reduces  to  a  constrained  root¬ 
finding  problem  that  can  be  solved  using  clipping  logic 
as  a  pre-processing  step  coupled  with  standard  root¬ 
finding  techniques.  Alternatively  nonlinear  program¬ 
ming  techniques  can  be  used,  although  it  is  likely  that 


Figure  3:  Clipping  of  Infeasible  Moment  Command 

nonlinear  programming  may  be  too  slow  for  on-line 
control  allocation.  In  order  to  use  root-finding  tech¬ 
niques,  one  must  first  determine  if  a  moment  command 
is  feasible  or  lies  within  the  AMS  boundary.  If  the 
command  is  feasible  standard  root-finding  techniques 
can  be  used.  If  a  command  is  not  feasible,  it  must  be 
clipped  at  the  AMS  boundary  and  the  axis  prioritiza¬ 
tion  method  presented  in  the  previous  section  provides 
a  means  for  performing  this  operation. 


4  AMS  Boundaries  for  Multiple  Effectors 

The  generation  of  an  AMS  for  a  single  LR  pair  is  rela¬ 
tively  straightforward.  Here  we  would  like  to  develop  a 
method  for  generating  the  AMS  for  multiple  LR  effec¬ 
tor  pairs.  To  illustrate  how  this  might  be  accomplished 
we  consider  the  rolling  and  yawing  moments  that  can 
be  generated  by  two  LR  effector  pairs  of  the  form 

Nt  =  Ni6 1  N2S2  + 

Lt  =  +  L2S2  +  ^3(53  +  ^4^4, 

and  for  the  purposes  of  example  we  let  Ni  —  —N2  = 
34/t  —  Ih/deg^^  Ns  =  —N4  =  iVi/25,  Li  =  — L2  = 
-4610ft  -  Ib/deg,  and  L3  =  -I/4  =  Li/5.  We  set  all 
effector  deflection  limits  to  ±30^6^.  To  generate  the 
upper  AMS  boundaries  for  the  1-2  effector  pair  we  set 
-  S4  =  0,  Si  =  30  =  -30,  and  let  62  vary  con- 

tinously  between  ±30.  The  lower  boundaries  and  the 
boundaries  for  the  3-4  effector  pair  can  be  generated 
in  a  similar  fashion.  To  generate  the  AMS  bound¬ 
ary  for  both  of  the  LR  pairs  acting  together,  one  may 
add  the  AMS  of  the  3-4  pair  to  each  point  on  the  1-2 
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AMS  boundary.  Figure  4  shows  the  effect  of  perform¬ 
ing  this  operation.  In  principle,  the  AMS  boundary 
could  be  generated  by  exhaustion  by  moving  the  ori¬ 
gin  of  the  3-4  AMS  boundary  to  each  point  on  the  1-2 
AMS  boundary  and  drawing  a  composite  AMS  for  the 
1-2-3-4  effectors  by  drawing  a  curve  that  touches  the 
extremal  points  in  the  roll-yaw  plane.  Such  a  proce¬ 
dure  would  clearly  be  inefficient. 


Figure  4:  Generation  of  AMS  by  Superposition 


Figure  5:  Generation  of  the  AMS  Boundary 


the  composite  boundary  as  those  that  produce  maxi¬ 
mum  or  minimum  values  of  the  h{x,x)  as  x  is  varied 
over  its  range  of  physically  realizable  values. 

Definition  2  The  composite  boundary  b{x)  gener¬ 
ated  from  all  physically  realizable  values  of  h{x,x)  = 
g{x  —  x)  +  f{x)  is  defined  as  the  set  of  extremal 
values  of  h{x,x)  for  {x  G  ^  <  ^max/} 

^X  G  "b  ^mirig  ^  ~  ^  ^max/  “b  ^maxg  } 

where  Xfninf^^ming  o<nd  Xfnaxfi^maxg  mini¬ 

mum  and  maximum  values  of  x  and  x  for  which  f{x) 
and  g{x)  physically  exist. 


It  is  therefore  natural  to  ask  if  there  exists  a  method 
or  condition  that  relates  the  AMS  of  each  LR  pair  con¬ 
sidered  individually  to  points  on  the  composite  bound¬ 
ary  approximated  by  the  outline  of  Figure  4.  The 
condition  for  a  point  on  the  boundary  turns  out  to 
be  related  to  the  slope  of  the  AMS  boundary  curves. 
For  convenience  we  shall  consider  the  top  left  quad¬ 
rants  of  the  two  boundaries  and  denote  the  bound¬ 
ary  curve  for  the  first  LR  pair  f{x)  and  the  bound¬ 
ary  curve  for  the  second  LR  pair  g{x)  as  shown  in 
Figure  4.  In  order  to  generate  the  composite  bound¬ 
ary,  one  can  generate  a  number  of  auxiliary  functions 
h{x,x)  =  9{x-x)-\-f{x)  where  the  point  {x,  f{x))  is  a 
point  at  which  the  boundary  curve  for  the  second  LR 
pair  attaches  to  the  boundary  curve  for  the  first  LR 
pair. 

Definition  1  An  attachment  point  to  f{x)  is  a  point 
P(x,  f{x))  to  which  the  origin  of  a  curve  g{x)  is  trans¬ 
lated  for  the  purpose  of  calculating  values  of  a  new 
function  that  is  composed  of  the  two  original  functions 
h{x,x)  =  g{x-x)  -f-  f(x) 


Theorem  1  states  a  necessary  condition  for  a  point  to 
lie  on  the  composite  boundary. 


Theorem  1  If  f{x)  and  g{x)  are  with  monotonic 
slopes  and  h{x,x)  =  g{x  -  x)  +  f{x)  then  for  any  x  G 
{x  G  H-  XmirXg  ^  {x  -  a;*)  <  Xmaxf  "b  Xffiaxg^ 

a  point  on  the  composite  boundary^P{x,h{x^x*))  sat¬ 
isfies  the  following  condition 


d{x  -  x) 


^9{x  -  x) 


X=X* 


(5) 


Proof:  By  definition  points  on  the  composite  bound¬ 
ary  are  extremal  values  of  h{x,x)  =  g{x~-x)-\-  f{x).  It 
is  assumed  that  f{x)  and  g{x)  have  continuous  deriva¬ 
tives  and  that  these  derivatives  are  monotonic  func¬ 
tions  of  X.  At  any  given  value  of  x,  we  are  interested 
in  finding  an  attachment  point  P(x*,/(x*)  that  will 
result  in  an  extremal  value  for  the  composite  function 
h{xjx).  The  first  order  necessary  condition  for  an  ex¬ 
tremum  of  h{x^x)  is 


Hence,  h{x,x)  shifts  the  origin  of  the  original  curve 
g{x)  to  a  point  on  f{x)  at  ar  =  x.  We  define  points  on 


=  0. 
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lx=x* 


(6) 


Applying  the  chain  rule  for  derivatives  and  after  rou¬ 
tine  manipulations,  it  is  easy  to  verify  that  Equation  6 
is  equivalent  to  Equation  5.  Since  the  slopes  of  f{x) 
and  g{x)  are  monotonic,  there  is  one  and  only  one 
value  of  X  that  satisfies  Equation  5  and  this  value  of  x 
is  defined  asx*,  ■ 


J^eixiark  1  If  (^x  3?  )  ^  ^maxg  thcji  the  bowfidavy 

point  is  given  by  P[x*  +  Xmaxg,h{xmaxg,'^*)]- 


Remark  2  If  (x  - <  Xmirxg  then  the  boundary 
point  is  given  by  P\x*  +  Xmirig  >  h{XmingiX  )]. 

Theorem  1  and  Remarks  1  and  2  collectively  state  that 
the  point  on  h{x,x)  that  has  the  same  slope  as  f{x) 
at  the  attachment  point  is  a  point  on  the  composite 
boundary,  as  long  as  that  point  does  not  require  vio¬ 
lation  of  the  effector  position  limits.  Theorem  1  may 
be  used  as  part  of  a  process  to  determine  the  com¬ 
posite  boundary  by  solving  the  Equation  5  for  x*  as  a 
function  of  x  only.  Contributions  to  the  AMS  from  ad¬ 
ditional  LR  pairs  can  be  accommodated  by  repeatedly 
applying  this  procedure.  A  similar  procedure  should 
provide  a  means  for  computing  AMS  volumes  in  mo¬ 
ment  space;  however,  the  planar  case  was  explored  here 
for  concept  development. 


5  Composite  AMS  Calculation 


This  section  provides  an  example  of  how  to  calculate  a 
section  of  a  composite  AMS  boundary  using  Theorem 
1.  It  is  relatively  easy  to  to  generate  segments  of  the 
AMS  boundary  for  a  single  pair  of  left-right  effectors. 
Theorem  1  provides  a  means  for  computing  composite 
boundary  segments  given  the  AMS  boundaries  for  sin¬ 
gle  L-R  pairs.  We  will  now  consider  the  case  where  the 
rolling  and  yawing  moments  produced  by  two  sets  of 
LR  pairs  as  in  Equation  4.  Note  that  even  though  the 
contributions  of  the  individual  effectors  to  the  yawing 
moment  are  non-linear,  there  is  no  coupling  between 
the  effectors.  This  implies  that  the  total  rolling  and 
yawing  moments  that  are  achievable  by  using  the  1- 
2  and  3-4  effector  pairs  together  can  be  obtained  by 
summing  the  contributions  of  the  respective  LR  pairs 

Nt  =  Ni2+Ns4, 

Lt  =  -Li2  +  L34. 


In  this  example  we  consider  the  following  parametric 
equations  for  the  segment  of  the  AMS  boundary  for 
the  1-2  effector  pair  alone  in  the  second  quadrant  of 
the  roll-yaw  plane 


-^12  —  ’ 
Li2  =  Kl8l  -  Kl62rr,a:.  > 


where  Kn  =  — 34/f  —  Ibldeg*^,  Ki  =  4610/^  —  Ihjdeg 
and  52m  ax  =  ^ote  that  all  effector  position  limits 
are  ±30®  in  this  example.  Equation  8  will  generate  a 
curve  of  constant  effector  deflection  for  52maa  = 
when  Ni2  is  plotted  as  a  function  of  L12.  This  curve 
is  a  segment  of  the  AMS  boundary  for  the  1-2  effector 
pair.  Eliminating  the  parameter  Si,  and  solving  for 
A^i2(-^'12)  one  obtains 

iVl2(l>12)  =  (9) 

Similarly  one  can  generate  an  equation  for  the  AMS 
boundary  for  the  3-4  effector  pair  acting  alone.  Again 
considering  the  AMS  boundary  in  the  second  quadrant 

iV34(L34)  =  (10) 

where  (7„  =  iifn/4  and  Q  =  Ki/2.  Equations  9  and  10 
can  be  plotted  in  the  roll-yaw  plane  and  for  conve¬ 
nience;  we  eliminate  the  subscripts  on  L  and  treat  it 
as  the  independent  variable.  Physical  limits  on  L  must 
be  observed  for  the  1-2  and  3-4  pairs  in  order  to  test 
the  conditions  stated  in  Remarks  1  and  2.  These  limits 
are  given  by 


^12ma.  - 

i^l2m.n  =  Wlmin  =  -276600/^  -  lb, 

=0, 

=  ^CiSlmin  =  -138300/i  -  lb. 

We  now  apply  Theorem  1  and  Remarks  1  and  2.  First 
we  arbitrarily  choose  the  1-2  cu^e  as  the  base  curve, 
which  we  will  refer  to  as  Ni2{L),  The  3-4  bound¬ 
ary  curve  Nz4  (L)  is  chosen  as  the  attaching  curve  and 
the  origin  will  in  principle  be  shifted  to  all  points  on 
Ni2(L)  viz.  iV34(L““L).  Theorem  1  states  that  points 
on  the  composite  boundary  occur  where  the  slopes  of 
the  base  and  attached  curves  are  equal.  The  slopes  of 
Equations  9  and  10  are  given  by 

=  2^(1  +  K,S2,„^J,  (12) 

^7{L-W  " 

Setting  Equations  12  and  13  equal  and  solvingjor  the 
value  of  L  that  satisfies  the  equality,  namely  L  ,  one 
obtains 


L  —  — — 


Kn  I  Qj]l 


(14) 


Now  that  L*  can  be  written  as  a  function  of  L  only, 
one  can  write  the  equation  for  the  composite  boundary 


(8) 
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(without  regard  to  position-limit  imposed  constraints 
on  L  -  L*)  as 

Nb{L,L  )  =  ^  -\-2-^S2^^^L 

+^(L-r)2+2^52_(L-r). 

(15) 


Remarks  1  and  2  must  now  be  applied  to  ensure  that 
it  is  physically  possible  to  move  the  control  surfaces  to 
points  where  the  slopes  are  equal.  If  this  is  not  phys¬ 
ically  possible,  due  to  position  limits  on  the  effectors, 
the  boundary  point  will  be  evaluated  at  the  upper  or 
lower  limit  of  L  for  the  effector-pair  in  question.  In 
this  case  the  attaching  curve  if  the  3-4  AMS  bound¬ 
ary  whose  upper  and  lower  bounds  on  L  are  given  in 
Equation  11.  Applying  Remark  1  we  find  that  when 
L  -  IT  >  5  bouudary  becomes  a  function  of 

L  only 


Kn-r^ 


Nb{L  )  =  j^L  + 


2  ■_  2— Jo 


Ki 


(16) 


Similarly  from  Remark  2  we  find  that  when  L  —  L  < 
LzAmin 


NBir)  =  ^L*^  +  2^S2„ 


(17) 


Figure  5  shows  the  composite  boundary  segment 
Nb{L^T,*)  for  the  example  above.  Construction  of  the 
other  boundary  segments  will  follow  a  similar  proce¬ 
dure  and  opportunities  to  exploit  symmetry  obviously 
exist  but  were  not  considered  here  for  brevity.  Also, 
it  is  clear  that  as  additional  boundary  segments  are 
computed,  there  may  be  overlap  in  some  regions  and 
a  method  must  be  devised  to  choose  the  most  liberal 
boundary  segment  at  a  given  point. 


6  Conclusion 

Concepts  for  new  constrained  control  allocation  strate¬ 
gies  were  explored  in  this  paper  that  deal  with  systems 
where  moments  are  nonlinearly  related  to  effector  de¬ 
flections,  such  as  those  encountered  in  the  case  of  yaw¬ 
ing  moment  contributions  from  left-right  effector  pairs 
on  aircraft.  The  concepts  were  illustrated  by  consider¬ 
ing  single  and  multiple  left-right  pair  effector  mixing 
problems  for  moments  that  lie  in  the  roll-yaw  moment 
plane.  The  attainable  moment  sets  for  individual  left- 
right  effector  pairs  are  relatively  easy  to  construct,  and 


Figure  6:  Composite  Boundary  Segment 


a  method  was  presented  that  allows  one  to  construct 
a  more  complex  composite  boundary  from  the  individ¬ 
ual  left-right  pair  attainable  moment  sets.  Necessary 
conditions  for  points  on  the  boundary  of  the  attain¬ 
able  moment  set  for  multiple  left-right  effector  pairs 
were  derived  for  the  planar  case  where  only  two  mo¬ 
ments  are  specified.  Extension  of  the  theory  to  include 
necessary  conditions  for  points  on  the  boundary  of  the 
attainable  moment  set  in  3-dimensional  moment  space 
presently  remains  an  open  but  important  problem.  A 
method  for  clipping  unattainable  moment  commands 
with  axis  prioritization  was  presented. 
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